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EFFECT OF RANDOM FLUCTUATIONS IN PRESSURE 
GRADIENT ON CBANNEL FLOW 
Ly Morris Perlmutter 
Iiewls Research Center 

National Aeronautics and Space Administration 
Cleveland^ Ohio 

ABSTRACT 

A randomly fluctuating pressure gradient of a stationary Gaussian Markovian 
form will cause a randomly fluctuating velocity to he superimposed on the steady 
Incompressihle flow in a channel. Correlations^ spectra^ and frequency re- 
sponse functions for the random functions are given. Random pressure signals 
were generated using Fourier series expansion with coefficients randomly 
picked from distributions whose parameters were obtained from the spectra of 
pressure signal. The random velocity signals were then obtained from the 
pressure signal by use of the frequency response function calculated from the' 
equation of motion. The increased power loss due to the fluctuations is 
given and the random pressure and velocity signals are compared for amplitude^ 
frequency^ and time lag. 

INTRODUCTION 

In many practical situations, flows in channels will have random pressure 
gradient fluctuations. It is of interest to understand the effect of the ran- 
domly fluctuating pressure gradient on the characteristics of the flow. The 
analysis is restricted to channel regions sufficiently downstream from the 
channel entrance so that entranee effects can be neglected. The analysis 
gives the mean square value of the velocity fluctuation across the channel 
as well as the increased power loss in pumping the fluid due to the fluctu- 


ations in the velocity. 
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A random sampling or Monte Carlo approach is also nsed^ :to generate a 
specific pressure gradient signal having a given power spectrimi and a 
Gaussian distribution of amplitudes. This calculation can he rapidly carried 
out hy use of fast Fourier transforms. The randomly varying pressure gra- 
dient signal generated is then used to find the corresponding fluctuating 
velocity signal at various positions across the channel. These nmerical 
signal: results can he useful^ for instance^ in numerical convective heat- 
transfer calculations in which numerical values for the randomly fluctuating 
velocity are required as input to calculate wall heat transfer. 


SBIBOL 


complex coefficients 


spacing between parallel plates 
frequency increment, ^/^p 
maxlmixm frequency 
frequency response function 
gain factor 


weighting function or impulse response function 


number of samples 
pressure gradient 
correlation (x, y) 


mean square power spectral density 


period of record 
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t dimensionless time^ vT/d^ 

At time between sample 

U fluid velocity 

Ug steady dimensionless velocity, -vpUg/d^Pg 

V transient dimensionless velocity, -vpU^/d^Pg 

V coordinate across channel 

Yj^ 4 sin Ky/K 

y dimensionless coordinate across channel, Y/d 

y dimensionless pressure gradient, P-t/Pg 

A measure of rate of fluctuation of pressure signal 

1 dimensionless meassure of rate of fluctuation of pressure signal, Ad^/v 

V kinematic viscosity 

p density 

variance measure of pressure fluctuations 
T dimensionless time difference, tg - t^ 

cp phase factor 

CO angular frequency, 2jrf 

( ) ensemble average 

* complex conjugate 

Subscripts; 

F fundamental frequency or lowest frequency, l/Tp 

H related to frequency response function 

I imaginary 

k finite spatial sine transformed variable 
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R real 

s steady component 

t fluctuating component 

T related to pressure signal 

(JD Fourier time transformed varialile 

ANALYSIS 


For viscous incompressible flow "between parallel plates with constant 
properties in the fully developed flow region (see fig. l) the momentum equa- 
tion can be written as fref. l) 

5u/ST = - l/p P(T) + V 

We can let the velocity and the pressure gradient consist of a steady part 

and a nonsteady part 

U = U + U, : P = P + P , 
s t^ St 

Normalizing the variables as follows 


- vpUg/d^P V = - vpU^./d^P y = Y/d t = vT/d^ T = P^/Pg 

( 2 ) 

we obtain for equation (l) 

d\ /dy^ = 1 or u_ = y - y^/2^ Sv/St = t + d^v/5y^ (5a;b) 

Taking Fourier transforms of equation (6) we obtain 


V^lCD 




CO 


+ a yjhj 


( 4 ) 
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'Taking the finite sine transform of equation (4) gives 


= ly3^5. . . 


( 5 ) 


where K = kjt. By solving for v^ and taking the inverse finite Fourier 


transform^ we obtain 


'CD 


Ta^ 


( 6 ) 


where 




H 


O - icPrr 

Yj^(k 2 + ico)= Hr + iHj = lH|e » 


(7) 


k=l,3^5 

where Yj^ = 4 sin Ky/K. The H is sometimes called the frequency response 
function^ j H| the gain factor^ and cpg- the phase factor, where 


(HI=^ = + H^j cpjj = tan-l(Hi/HR) 


(8;9) 


The gain factor is plotted in figures 2 and 3; the phase factor is 
plotted in figure 2 and is discussed more fully later where used. We can 
solve equation (6) for v hy taking the inverse transform. Using the con- 
volution theorem gives 


V = 


I 


+00 /•+<» 
gi03t jj (if ^ 

-00 f -00 




h(9,y)r(t - e)de 


(10) 


where h is the weighting function given hy 

r oo 


He) = 


r He^- 

J-” 


Y^e 


-K^(e) 


for 9 > 0 


df = < 


k=l,3,5 


) ( 11 ) 


0 


V . 


for 6 > 0 


J 
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The h. and H are a Fourier transform pair. Since the process is station- 
ary^ we can write the cross correlation equation as 


(r(ti)v(t3^ + t)) = R.py.(T) 



h(0,y)RyY('f - 0)d0 


( 12 ) 


where " ^1 ” cross spectrum is the Fourier transform of the 

cross correlation, hy using the convolution theorem as before, we can write. 

= H(oj)S^(oo) (13) 

We can find the velocity autocorrelation as 





h(0,y)R^y(T 


+ 0)d0 


(14) 


We can then write the velocity power spectrum as 


S^(co) = IH(od)1^S^(oj) (15) 

This shows the velocity spectrum is equal to the pressure spectrum 
times the square of the gain factor. 

Pressure Gradient Correlations and Power Spectrum 
A common form of random fluctuations encountered in physical systems 
is called a stationary Gaussian Markoff random process. Markoff processes 
are random processes whose relation to the past does not extend beyond the 
immediately preceding observation. As is shown in reference 2, page 215, 
the autocorrelation function for this random process is given by 

R^ = 0^e = 0^e~^- where = R^( o) = (y^) (16) 

The rate of fluctuation A and the time span T have been nondimensionalized 


to be 


A = A d^/v ; T = Tv/d^ 
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The Fourier transform of Ryy gives the power spectrum of the pressure 
fluctuations as a simple Lorenzian form: 

= 2Aa^/(,A^ + (17) 

p 

The cr^ is a dimensionless measure of the amplitude of the random fluc- 
tuations (t^)- The larger the value of 0^^ the larger the amplitude of the 
pressure fluctuations. The A is a dimensionless measure of the average 
number of fluctuations per unit time. It can also he considered as^ the in- 
verse of the dimensionless characteristic decay time of the autocorrelation. 
The larger the value of A, the greater the fluctuation rate of the signal. 
The power of the pressure signal is given hy the integral of the power spec- 
trum over all frequencies. This can he written as 


% 


^+00 

Ryy-( O) = I S. 

J “00 


rr 


df = aj = <r^) 


(18) 


The normalized power spectrum for the pressure gradient fluctuations 
Syy/(r^) is plotted in figure 3 for different values of the rate of fluctu- 
ation parameter A. It can he seen that for smaller values of A much more 
of the spectral power is located in the lower frequencies since the frequency 
hand for A of 0.1 is much narrower than in the case for A = 100. 

Fluctuating Velocity Mean, Correlations, and Power Spectrums 
The pressure-velocity cross correlation for the present case can he 
written using equations (ll), (12), and (16) as 

r 




r 



[Yj^/(k2 - aJ e-^^ -[2A/(A + K^jj, 


2 

.-K T 


for T > 0 




k=l,3,5 


j ^ + a) for t < 0 

^ k=l,3,5 


( 19 ) 
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Also the cross spectrum can he written using equations ( 7) ^ (13)^ and (17) as 


=1 ^ Yj^/(k 2+ icn)\,2Aa^/(A2 + oD^) 


-k=l,3,5 


In a similar manner we can calculate the velocity autocorrelation from 


equations (ll)/il4), and (l9) as 


CO 00 




(2a/[a + L^])e“^^V(K2 + 


k=l,3 7=1,3 


Similarly, we can write the velocity fluctuation power spectrimi using equa- 
tions (8), (is), and (17): 

S^(iu) = I H(cu) 1 2 s^(cd) = (i^ + h|)s^ (22) 

From equation (lO) we can see that the ensemble average of the fluctu- 
ating velocity (v) at any given y is zero since the ensemble average of the 
fluctuating pressure gradient (t) is zero. 

The mean square velocity (v^) can be obtained from equation (21) since 
E^(o, y) = (v2). These results are shown in figure 4. Notice that the mean 
square velocity is equal to the total power of the velocity spectrtim 


R (o^y) 

YV' 


S df = 
w 


|H(cD)rs^(co)df 


The spectral power or the mean square velocity is given by the product -of the 
square of the gain factor times the pressure spectrum integrated overall fre- 
quency. Prom figure 2 we can see that the magnitude of the gain factor is 
smaller near the wall (y = O.l) compared to the centerline (y = 0.5), so 
that the mean square velocity is smaller near the wall as compared to the 
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centerline. The mean square velocity goes to zero at the wall in figure 4, 
as would he expected from the zero wall velocity boundary condition. The 
gain factor in figure 2 is greatest at the lower frequency and has a high- 
frequency cutoff. Since the spectrum for the pressure gradient (fig. 3)^ is 
also greater at the lower frequencies for low values of A as compared to 
the high A case^ the mean square velocity will he larger for the lower 
values of A, as can he seen in figure 4. 

The fluctuating velocity power spectrum given hy equation (22) normal- 
ized hy the total power of the spectrum at a given position in the channel 
Sy^(y_, f )/Rvv(y^o) is also shown in figure 3. These curves are discussed 


more fully later. 

Power Dissipation Due to Velocity Fluctuations 
It is desirable to know how much additional power must he supplied to 
the flow system to maintain the same average flow rate for the fluctuating 
flow as compared to the steady flow case. As shown in ref . (l), the ■ 
external rate of work done hy the pressure force w^ is equal to the sum of 
the rate of increase in kinetic energy K and the rate of dissipation of 
energy due to internal friction w^. If we then take the ensemble average^ 
the rate of change of the kinetic energy term will he zero because the 
process is stationary. We then have the external rate of work done hy the 
pressure force Wg equal to the rate of dissipation of energy due to ex- 
ternal friction w^ 


-(UP) = 



(24) 


To find the rate of dissipation of energy due to internal friction we 
can then write 

wj = -(UP) = -U3P3 - (U^P^) = Wf^g - (U^P^) 


(25) 
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If we let W|-^ g te the rate of dissipation due to steady flow^ the ratio of 
increased dissipation due to the unsteady flow over the steady flow dissipa- 
tion can he written as 

(Wf - a (U^P^)/tjgPg = 12(^) = 12R^^(o) (26) 

Integrating equation (l9) over y and substituting into equation (26) give 

CO 

1/K^(K^ + A) (27) 

k=l, 3, 5 

A plot of equation (27) is given in figure 5 as a function of the fluc- 
tuating rate parameter A. As A becomes very large^ that is^ fluctuations 
of the pressure gradient become very rapid, the frictional power loss re- 
duces to the steady power loss. The frictional power loss increases with 
smaller values of A or slower fluctuations of the pressure gradient. 

Generation of Random Signal by Model Sampling 

It would be helpful in certain cases to be able to generate a random 
signal analytically on the computer consistent with the power spectrum of the 
signal. This would be useful, for instance, in obtaining numerical solutions 
or various other statistical characteristics of linear or nonlinear equations 
of random systems that may be difficult to obtain by more usual analytical 
procedures. A method for generating random signals is given in peferen.ee 3. 

A somewhat different approach is used in the present analysis. 

Discrete Fourier transforms . - When a waveform is to be analyzed on a 
digital computer, it is the discrete Fourier transform that must be used. 

This can be derived from the usual Fourier transform as shown in reference 2, 


page 56, as follows. 
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If T(t) is periodic in T . we can write the Fourier expansion 

Sr 


where 


^ in2jtt/T 

r(t) = ^ a^e 


= 1/Tp f ^ dt 


¥e can see from equation (28) th^t 


lim(2?pan) - 

Tp-«o 


Similarly, if is periodic in f^, we can write 


00 


-i2Ttnf/f 


where 


dn ='¥fp)yj, ,2 

We can see that, similarly, 

lim(f dj'= T(t) (2 

fp-^ ^ 

If we take the time between sample values of T-^ as At, we can take 
N samples so that W At = Tp. Then, assuming that At = l/fp we can 
write equation (28) as 


±2nf/t 


j,^k At) = 2^ 

n=-lJl/2 
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¥e caxi also write equation (31), assiming that N Af = f , as 

N/2 

=(i/N) V (35) 

M/2 

Since can he related to the Fourier transform hy equation (so), 

equations (34) and (38) form a discrete Fourier transform pair that is an 
approximation of the Fourier integral transform pair. 

A mmierical method for evaluating either equation (34) or (35) in a ^ 
very efficient manner on digital computers has heen developed called the fast 
Fourier transform (ref. 4). One implementation of the fast Fourier trans- 
form is available as an IBM-scientific subroutine package called HARM/DHARM. 
Many others are also available. 

Gaussian random process . - Since Tjj. is a real function a^ must be 
an even function and j must be an odd function around n = 0. Then we 
can write 


^k = 


S'n R ‘^os(2whkyN)- j sin(2rtnk/w) 

n=-N/2 


Us) 


We can assume a Gaussian random process for This process is commonly 

used as an idealization of many natural phenomena associated with the super- 
position of a large number of many small effects. It is shown in refer- 
ence 5, page 160, that for a Gaussian random process the following condi- 
tions are true: a is a random variable with a Gaussian distribution also 


(®‘n,R^k,R) - ^®'n,I®^k,l) = 



k=n 


(anbk) = 0 


(37;38) 


0 k^n 

Then, spectrum of T(t) can be shown to be related to Un) 
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Ia>=n2jt/T 


= 2T. 


when 


n f 0 


(39a) 


when 


n = 0 


(39h) 


Model sampling to generate signal . - To generate a random signal T(t) 
we can use equation (34) and the fast Fourier transform program. The terms 
and an^I are random variables randomly picked from the joint Gassian 


distribution 


l) = (1/2" (4 ) 


Following the method given in reference 6, page 39, we can randomly choose 
values of a^ p and I "tliat will satisfy the distribution of equa- 
tion (40) by using the following, equations: 


%R = (2(9^)) (-In cos 2nRg 


where Rg and Rj. are two different random numbers randomly picked from a 
uniform distribution between 0 and 1. The values of R can readily be gen- 
erated by the computer. The values of (a^j can be obtained from equa- 
tion (39) using the appropriate power spectrum S^y. Recalling that 
is an even function and I is an odd function, in the particular fast 
Fourier program available (HARM/dHARM), we only need find a^ and a^^ j 
from n = 0 to N/2. Then a(R/2)+l R “ ^(n/ 2)-1 R^ Similarly, for 

a^ j values were only needed for n = 0 to n = N/2; since S'(]\r/2)+l I 
sc^iJialod. 1 I-^ cbc*^ s<Q -j- 


was equal to zero- 
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Pressure Gradient Signal Results 

The normalized fluctuating pressure gradient T is shown in 

figure 6 for different values of A. The pressure gradient power spectrum 

used in obtaining these results was assumed to he a stationary Gaussian 

Markoff process and is given hy equation i[l7) and plotted in figure 3. The 

results were calculated for values of N of 128 and AT taken as 100. 

P 

The curves were computer plotted hy joining the output points hy dotted 
lines. Since the spectrum for the higher fluctuating rate parameter A = 100 
is greater at the higher frequencies^ as can he seen in figure '5, the re- 
sulting pressure signal for large A has many more fluctuations and cross- 
ing of the 0 value line per unit time than is the case for the lower values 
of A. 

Generating Velocity Signal From Pressure Signal 
The method of calculating the velocity signal from the pressure signal 
is as follows. We can write, similar to equation (36), 

00 

cos 2jtnk/n + hj^ j sin 27tnk/w) (42) 

n=-co 

where h^^ is an even function and h^.^ j is an odd function. From equa- 
tions (6) and (30) we can write 

~ "'^n'^p ~ (^jo:)=n2jr/Tp|^®’n^p) 
or 

h^ = 

Thus, hy finding the values of a^ for a given signal discussed 

previously, we can use equation (43) to obtain the new values hj^ which can 
he used in equation (42) to obtain the the -velocity signal - v :using'.the fast - 
Fourier transform program. 
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¥e can see from equation ( 4S) that the mean values of and I 

are zero. Also^ since and aj^^I are normally distributed^ ^n^R 

hn I are also normally distributed. The variance of b^^ can be seen to be 

("^n.R) = + ^,l) (^n)l®uP 

except at n = 0 

{<^)n,0 - ° 

('>o,h) > («o,r)Ho,E 

The spectrum in this case would be given as in equation (39) by 

S„(n) = Tp2(b^) = Tp2(4)|Hl^ = V")Kr 

except at n = 0 when 




This shows that b^^R and b^^I cannot only be found from a^ hut 
can be picked directly from a normal distribution with a variance of b§ 
given by equation (46) as in the case of a^^ values. 

In figure 6 the normalized fluctuating velocity at the centerline of 


the channel Vy=o. 5 is plotted as solid lines. These results 

were calculated from equation (42) using the fast Fourier transfom routine. 
The coefficients b^ were found from equation (43)^ which gives the values in 
terms of an which are the coefficients of the random pressure fluctuations. 
These were randomly chosen by equation (4l). Thus, the velocity shown by 
the solid line on each graph is caused by the pressure fluctuation shown as 
the dotted line on the same graph. Notice that all the curves have been 


normalized by their own root mean square values. 
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For low values of Tk, the shapes of the curves are very similar hut 
there is a small time lag in the velocity curve compared to the pressure 
curve. This time lag can he explained hy postulating a distortionless fre- 
quency response function H' with a distortionless gain factor [H(cja)| = |H| ' 
and a distortionless phase factor cp' = -oitQ. 

Then equation (6) can he written as 


, , , icut„ 


r. 


CD 


Taking the inverse transform and using the time shifting theorem gives 
v(t) = |H[ ’r(t + to) 


X48) 


(49) 


This shows, for the distortionless response function, the velocity signal 
lags the pressure signal hy a time increment to^ while the amplitude is 
changed hy an amount |H|-. 

In comparing the distortionless frequency response function to the 
present response function we can see from figure 2 that for small angular 
frequencies the response function approaches the distortionless case. We 
can calculate the time lag to for small angular frequencies hy solving for 
cp as CD -»• 0. This gives, from equation (9), cp^^ = -cd(0. 0781). Thus, the 
time lag between the velocity and pressure signals for small cd is given hy 
Iq = 0.0781. Thus, in figure 6 for A = 1 we see the velocity signal is 
close to the pressure signal with a small time lag between them. This is 
because the response function is close to distortionless for small cd and 
the pressure signal for small A is concentrated in the small cd region 
(see fig. 3) . 

For the higher values of A we can see that the velocity signal does 
not contain many of the higher frequency fluctuations that are present in the 
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velocity signal. This is because the gain factor (see fig. 3) has a high- 
frequency cutoff which caused the velocity spectrum to have a similar high- 
frequency cutoff. Since at high A much of the pressure spectrum is in the 
higher frequencies^ these higher frequencies are not transmitted to the ve- 
locity signal. 

In figure 1 , the velocity signal at a given position y in the channel^ 
normalized by its root mean square value at the same y value^ is plotted at 
two different values of y for the same pressure signal. For low values of 
A the curves near the center of the channel (y = 0.5) are very similar to 
that nearer the wall (y = 0. 01) . This happens because the transfer function 
for small cn is close to distortionless. The velocity signals near the wall^ 
although of a smaller amplitude than the velocity signal near the center of 
the channel (see normalizing factor in fig. 4)^ both have almost identical 
shapes. At the higher values of A we see that the centerline signal fluc- 
tuates much less than the velocity signal near the wall. This is due to the 
fact that the gain factors for y close to the wall are flatter and do not 

V- 

have as strong a high-frequency cutoff as the gain factors near the center. 
This causes the velocity spectrimi near the wall (y = O.l) given in figure 3 
to be flatter and extend to higher frequencies than the velocity spectrum 
near the center of the channel. 

At the higher values of 1 \, we can also notice a time lag between the 
signals. This effect is due to the phase factor shown in figure 2. It can 
be seen that the phase factor curves for the signals closer to the wall have 
smaller slopes and so have a smaller delay time than the velocity near the 


center of the channel. 
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RESULTS AMD CONCLUSIONS 

The effect of superimposing a randomly fluctuating pressure gradient of 
a stationary Gaussian Markovian form with a zero mean value upon a steady 
pressure gradient will cause a randomly fluctuating velocity component to he 
superimposed over the steady flow in a channel. The mean value of the ve- 
locity fluctuation is zero. The mean square value of the velocity fluctu- 
ations which is a measure of the amplitude of the fluctuations will he high- 
est in the center of the channel and reduce to zero at the wall. Also^ the 
slower the rate of the pressure fluctuations^ the greater the amplitude of 
the velocity fluctuations. 

The frictional power loss in pumping the fluid is increased hy the pres- 
sure fluctuations^ the slower pressure fluctuations giving larger power losses. 

The normalized velocity signal values were very similar to the normalized 
pressure signal hut with a dimensionless time lag tg of about 0. 0781 for low 
pressure fluctuation rates. At the higher pressure fluctuation rates the 
higher frequencies of the pressure signal did not appear in the velocity sig- 
nal. 

Comparison of the velocity signals across the channel normalized hy its 
local root mean square value showed that for low fluctuation rates in the 
pressure signal the values of the normalized velocity near the wall were 
almost Identical to the normalized velocities near the center of the channel. 
However^ at the higher pressure fluctuation rates the higher frequency fluctu- 
ations present in the velocities near the wall do not appear in the velocities 
near the center of the channel^ and the velocities near the center of the 
channel lag the velocities near the wall. 

The random pressure signal function was generated hy a Fourier series 
expansion where the coefficients were randomly chosen hy model sampling from 
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a frequency distribution whose parameters are given by the power spectrum of 
the signal. By use of the fast Fourier transform method of computation the 
function can be rapidly evaluated from its Fourier coefficients. The Fourier 
expansion of the fluctuating pressure signal can then be used to calculate the 
randomly varying velocity at various positions from the wall. These velocity 
results can be useful, when, the details. '.of a randomly fluctuatihg velocity 
distribution in a channel are needed^ as for instance in convective heat- 
transfer problems. 
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Figure 1. - Parallel plate channel flow model. 
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Figure 2. - Frequency response function of system. 
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Figure 4. - Normalized mean square fiuriuating velocity com- 
ponent across channel. 
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Figure 5. - Percentage Increase in power dissipation ratio due to 
flow fluctuations. 
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